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Abstract 



The phenomenon of six degrees of separation is an old but interesting problem. The con- 
siderations of the clustering coefficient reflecting triangular structures and its extension to 
square one to six degrees of separation have been made[T]. Recently, Aoyama[2] has given 
' some considerations to this problem in networks without loops, using a sort of general formal- 

ism, "string formalism". In this article, we describe relations between the string formulation 
proposed by Aoyama and an adjacent matrix. Thus we provided a reformulation of the string 
^ ' formulation proposed by [2] to analyze networks. According to it, we introduced a series of 

C/3 , generalized q-th clustering coefBcients. The available rules between diagrams of graphs and 

jy^ ■ formulae are also given based on the formulation. Next we apply the formulation to some sub- 

O ' jects in order to mainly check consistency with former studies. By evaluating the clustering 

* ^ , coefficient for typical networks studied well earlier, we confirm a validity of our formulation. 

Lastly we applied it to the subject of two degrees of separation. 

^ . 

P> I keywords: Six Degrees of Separation, String, Clustering CoefRcient, Degree Distribution, 

Generalized Clustering Coefficient 

■ 1 Introduction 

o ■ 

00 ■ In 1967, Milgram has made a great impact on the world by advocating the concept "six degrees 

I of separation" by a celebrated paper [3] written based on an social experiment. Six degrees of sep- 

■ aration indicates that people have a narrow circle of acquaintances. A series of social experiments 
made by him and his joint researcher [5] suggested that all people in USA are connected through 
about 6 intermediate acquaintances. Their studies were strongly inspired by Pool and Kochen's 
study [6J. At the time, however, numerically detailed studies could not be made because computers 
and important concepts, such as the clustering coefficient needed for a network analysis nowadays, 
have not yet developed sufficiently. 

One of the most refined models of six degrees of separation was formulated in work of Watts and 
' Strogatz0,[8]. Their framework provided compelling evidence that the small- world phenomenon 

is pervasive in a range of networks arising in nature and technology, and a fundamental ingredient 
in the evolution of the World Wide Web. Another is the scale free networks proposed by Barabasi 
et al.[2], [IQJ. Many empirical networks are characteristic future of scale free [H] , [13] , [II] , [TS] . In 
spite of furthermore exploring of six degrees of separation[TS|,[I7], they do not examine closely 
Milgram's original findings by their model, especially how influence can the clustering coefficient 
proposed in the paper ^ has. We have made some study of it in our previous paper [18^ by imposing 
a homogeneous hypothesis on networks. As a result, we found that the clustering coefficient has 
not any decisive effects on the propagation of information on a network and then information 
easily spread to a lot of people even in networks with relatively large clustering coefficient under 
the hypothesis; a person only needs dozens of friends. Moreover we devoted deep study to the six 
degrees of separation based on some models proposed by Pool and Kochen 6J by using a computer, 
numerically [15]. In the article, we estimated the clustering coefficient along the method developed 
by us jl8| and improved our analysis of the subject through marrying Pool and Kochen's models 
to our method introduced in JB]- As a result, it seems to be difficult that six degrees of separation 
is realized in the models proposed by Pool and Kochen[6j on the whole. 
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If the network of human relations has a tree structure without loops, a person connects new 
persons in power of average degree, when he(she) follows his(her) acquaintances step by step on 
his (her) network of human relations. Then the phenomenon of six degrees of separation is not 
so amazing, if a person has a few hundred acquaintances. A question is that networks of general 
human relations include some loop structures. This structures decrease the number of new persons 
that connected with him(her) when he(she) follow his(her) acquaintances step by step. One of 
indices characterizing loop structures is the clustering coefficient. Thus it will be important to 
investigate the effect of the clustering coefficient and degree distribution on the six degrees of 
separation. It is, however, difficult to investigate the influence of general loop structures. There 
are only a little research focused on the effect of loop structures. 

Newman first studied the influence of loop structures in a network on the subject jlj. The 
study is so stimulating but only triangle structures and quadrilateral structures on networks were 
considered. It seems to be difficult to generalize his framework to q-polygon. Recently Aoyama 
proposed the string formulation on the subject[2]. The idea inspired our study in this article, 
greatly. Unfortunately he considered only tree approximation in the structure of networks. Since 
he deals with mainly scale free networks, the approximation is valid up to a point. 

In this article we pursue the string formulation and try to discuss the infiuence of general loops 
to six degrees of separation. One of the aims in this article is to reformulate the string formulation 
based on an adjacent matrix. We can systematically analyze general networks with arbitrary loop 
structures by the reformulation. Next is to check the results derived from it are consistent with 
results studied so far. Then we apply it to the problem of two degrees of separation as a first step 
so that future prospect for the problem are opend. 

The plan of this article is as follows. After introduction, we reformulate the string formulation 
by an adjacent matrix in the following section 2. First we explain some notations used in this article 
and reformulate the string formulation by using an adjacent matrix. Here we introduce R matrix 
that play central role in our formalism. According the formalism, we introduce generalized q-th 
clustering coefficients as well as the usual global one. Next we give a diagrammatic interpretation 
for every term appearing in the expansion of the power of R matrix like Feynman diagram [21j 
in quantum field theory. In the next section 3, we evaluate some clustering coefficients on some 
typical networks and discuss the consistency with results investigated by now. We shall discuss 
the phenomenon of six degrees of separation in the section 4. Since any reliable conditions have 
not been given for p-th degrees of separation in networks with strongly connected components 
yet, we can not provide full discussions. Thus we discuss two degrees of separation and compare 
the results to those given by Aoyamap] in scale free networks. We find that the result is a little 
different from Aoyama's one [2]. The analysis of an adjacent matrix is seemed to support rather 
our results. We shall devote discussion on six degree of separation based on the formulation in a 
subsequent article [501. The last section 5 is devoted to summary. 

2 String Formulation and Adjacent Matrix 

We, basically, follow the string formulation introduced by Aoyama [5]. 
2.1 Notations 

In this section we describe notations used in this article. We consider a string-like part of a graph 
with connected j vertices and call it " j-string" . N is the number of vertices in a considering network 
and Sj is the number of j-string in the network. (Note that Sj in this article is N times larger 
than defined by Aoyama 2 .) By definition. Si = N and ^2 is the number of edges in 

the network. Sj is the number of nondegenerate j-string where a nondegenerate string is defined 
as strings without any multiple edges and/or any loops in the subgraphs as seen in Fig.l. We, 
however, define that the nondegenerate string contains strings homeomorphic to a circle. 

We call strings without any loops open string and strings homeomorphic to a circle closed 
string. Thus we consider closed strings and open strings in this article. 

It is so difficult to calculate Sj and Sj, generally. Aoyama has calculated Sj up to j = 7 but 
did not supply the explicit in his article j2j. His article says that it would needs dozens of papers if 
the full expression is described. It would be maybe impossible to calculate Sj and Sj with j > 7 
explicitly at the present moment. 
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[a) Degenerate strings 



[b) Nondegenerate strings 



Figure 1: Two types of strings 



2.2 Clustering CoefRcient 

By using the string formulation, we can defined the usual clustering coefficient which essentially 
counts the number of triangular structures in a network. Let be the number of polygons with 
q edges in a network. Following Newman [T], the usual global clustering coefficient C(3) is given by 



(3) 



6 X number of triangles 6A3 
number of connected triplets S3 



(1) 



We introduce some indices uncovering properties of general polygon structures except for triangle 
structures in a network. From the expression of Eq.(l), we can generalize it to q-th clustering 
coefficient C(^qj straightforwardly; 



(«) 



2p X number of polygons 2pAp 
number of connected p-plets Bp 



(2) 



2.3 Adjacent Matrix Formulation 

We reformulate C(g) introduced in Eq.(2) by utilizing an adjacent matrix A = {uij). Generally 
the powers, A^, A"^, A^, ■ ■ ■ of adjacent matrix A give information as to respecting that a vertex 
connects other vertices through 2, 3,4, • • ■ intermediation edges, respectively. The information of 
the connectivity between two vertices in A"' contains multiplicity of edges, generally. For resolving 
the degeneracy, we introduce a new se ries of matrices i?" which give information as to respecting 
that a vertex connects other vertices through n intermediation edges without multiplicity. We can 
find it by the following formula [22]: 



n (1 - 

p„ _ \ ^ ik,ij,ik-ij>l 



(3) 



where the product of (1 — Si^i^) of the numerator has the role of protecting of degeneracies from 
strings and (f — Si^i^) of the denominator is, however, needed to keep a closed string. 

This expression has (n — l)ply loops in a computer program and so it is almost impossible to 
calculate i?" within real time for large N when the rank of R is N. The expansion of Eq.(3) has 
2n(n-i)/2 |;gj-jj^g formally. This value is 32768 for n = 6 that is needed in the analysis of six degrees 
of separation as will be discussed in the later section. Though many terms really vanish, has 
still so complex expression. We only give the expressions of R^ ~ R^ here and will give that of R^ 
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in Appendix; 

[R%f^[A%f-{G,A},f+a.,f, 

[R%f = [A%f - {G.A^}^^ + {A,diag{A'')]^^ + 2\A\f + [G^ - G - AGA\; + ■ia,f\A\^ 
\R\i = [^']./ - {A,diag{A^))^^ {G, A^},/ - { A^, d*a5(A--')}^^ + 'i{\A%i)\A\,f 
+ 3[A3],/[A],^ + 2{G2, A},/ + \GAG\,s - 6{G, A},/ - {AGA,A},f + 3(^3]^^ 
+ {A,dmg(AGA)}^^ + 2{diag(yA^G%^ - {A ■ diag{A^) ■ A]if - [diag{A^)],f 

+ 3^a,kakf[[A%f + [A%k - S,f[A%f^ +4ai/(l - a,f), (4) 



where suffix is abbreviate in trivial cases and {•, •} means the anticommutation relation; {A, B} = 
AB — BA. diagA indicates the diagonal matrix whose elements are the diagonal elements of A, 
and G is the diagonal matrix defined by 



G = 



fci 
k2 
h 



(5) 



where ki is the degree of vertex i. 

By using i?", Sp and generalized p-th clustering coefficient G(p) are given by 



TtRp 



id 



-1' 



(6) 
(7) 



where denominator and numerator indicates the contribution from open strings and a closed string, 
respectively. Thus usual clustering coefficient G(3) becomes 



G, 



TtR-^ 



TrA3 



(3) 



\\A\\-TtA^' 



(8) 



where we introduced a new symbol || • • • || which denotes \\A\\ = J^i j Ai 



2.4 Diagrammatic Interpretation 

The expression of i?" is rapidly complicated as n increases. We, however, notice that every term 
appearing in the expansion of i?" closely corresponds to a certain graph p2]. So if a certain graph 
is given, we can write out the expression corresponding to it like Feynman's rule in quantum field 
theory We describe the diagrammatic construction of i?" based on the relation between both. 

1. Draw all graphs with n edges including degenerate graphs except for closed string. 

2. Assign sign (— for every graph where v is the number of vertices included in the 
graph. 

3. Calculate degeneracy index m defined in the following Eq.(9) and it is the coefficient of the 
term corresponding to the graph; 



m 



n 



hi 



1, 



\Ga, 



(9) 



where [ ]gq means Gauss symbol. The coefficient supply essentially the number of Euler 
paths in a graph. We do not distinguish a path and its opposite path and start from vertex 
with hight degree in the cases that there is not any odd vertex in the graph. 
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(a) (b) (c) 



Figure 2: Typical diagrains@I 




(a) (b) 



Figure 3: Typical diagrams II 

Thus the subject of finding the expression of i?" reduces to the estimation of diversity of graph 
with certain constant number of edges and the number of Euler paths in the graph. 

We describe the typical relations between expressions and graphs. The typical graphs arc given 
in Fig. 2 and Fig. 3 where the arrows indicate the order in Euler paths and i and f shows an initial 
vertex and a final one, respectively. These graphs correspond to the following expressions; 



Figure2{a) = h, (10) 

Figure2{b) = aif, (11) 

Figure2{c) = {A%i, (12) 

Figure3{a) = {A%f, (13) 

Figure3{b) = aij{A'^)if. (14) 



It is assumed that the right vertex has not any other edges in (a) of Fig 2. The generalization of 
Eq.(12) to polygons with n vertices is easily given by (^")ii. Eq.(13) is also straightly extended 
to {A'^)if in the cases that n vertices are linearly connected. 

Degenerate multigraphs also reduce multigraphs with little multiplicity as shown in Fig. 4. The 
left graph comes down to the corresponding expression Uif in (a) of Fig. 4. Since these are only 
useful relations between graphs and expressions, these graphs should be independently considered 
in the step 1 in the diagrammatic construction. These correspondence may be rather effective in 
estimating explicit expression of R^. 

2.5 An Example 

We give an example of the diagrammatic construction of R". As a nontrivial case, we consider the 
construction of R'^ which is constructed from graphs with four edges. All topologically independent 
graphs with four edges are given in Fig. 5. For every graph, signs and coefficients derived from step 
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(a) 

(b) 



Figure 4: Reducing of diagrams 




Figure 5: Topologically independent graphs with four edges 
2 and 3 in the diagrammatic construction are as foUows; 



Figure5{a) 




(15) 


Figure5{b) 


= [G%f, 


(16) 


Figure5{c) 


= -{A^,G}u, 


(17) 


Figure5{d) 




(18) 


Figureb{e) 


= -[AGA],,, 


(19) 


Figureb{f) 


= 2[A%f, 


(20) 


Figureh{g) 




(21) 


Figure5{h) 


= 3a^j[A%f. 


(22) 



It is easy to confirm that the sum from Eq.(15) to Eq.(22) give i?^ in Eq.(4). 

3 Clustering Coefficient of Some Networks 

In this section we calculate C(3) in some typical networks by using the formulation given in the 
previous section. Thus we investigate consistency of our formulation with results observed so far. 
We adopt the configuration model [23] . [24] , [25] in producing networks. The model can system- 
atically produce networks with arbitrary degree distribution. But the networks produced by the 
model are degenerate multigraphs, generally. We modify it a little to produce networks without 
multiplicity. Since it is not essential in this article, we omit the technical details of it. 

First we study random networks [26] where the degree distribution is Poisson distribution |5]. 
The N dependence on C(3) given by computer simulations is shown in Fig. 6. Theoretically it 
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♦ <k>=10 
■ <k>=25 
mi «k>=25) 



y = -1.0122X + 3.3049 
= 0.9998 



Log N 



Figure 6: The scaling of C(3) in random graphs with Poisson distribution 



♦ <K>=25,a=12 

■ <K>=25,a=5 

A <K>=10,a=5 

Mi «K>=25,g 

y = -1.0094X + 1.5729 
rf = 0.9984 




Log N 



Figure 7: The scahng of C(3) in graphs with Normal distribution 



is shown that the clustering coefficient behaves as C(3) = {k)/N in random networks where (fc) 
indicates average degree [H]. The log- log plot in Fig. 6 shows linear behavior with slope nearly 1. 
Thus N dependence of C(3) is consistent with observations so far. 

Moreover we can observe similar scaling law in networks with the normal distribution in degree 
distribution. As Fig. 7, all the slopes of the log-log plot indicating A^-C(3) relation are also nearly 
1 for various values of (fc) and the standard deviation a. The delicate difference of the values of 
the slopes was discussed in p2] . 

It is known that the clustering coefficient follows the scaling law C(3) ^ 7V^o.75 j^j. ^qqIq fj-gg 
networks[TT] with (k) = 4. Our simulation results for scale free networks that have the average 
degree 4 are given by Fig. 8. The slope, about -0.73, of the approximative line in Fig. 8 shows that 
the above result almost agrees with current one. This also strengthens the validity of presented 
formulation. 



5.3 



V — — O- 7 2 X — 0.4505 



F=t — O.S1 



. l_nC<:^3> ^mff& CI-nCCD3:)> 



Figure 8: The scaling of C(3) in scale free networks with (k) = 4 
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4 Application to Two Degrees of Separation 



Aoyama has proposed a condition, so-called Milgram Condition, for p-th degrees of separationj^]; 



M,^|^0(iV). 



For six degrees of separation, we obtain from Eq.(6) 



(23) 



(24) 



Before evaluating for six degrees of separation, we as a first step study two degrees of 
separation to check consistency in this article. Since 2-string cannot have any loops, we have only 
to consider graphs by the tree approximation in the study of two degrees of separation. 

Fig. 9 and Fig. 10 show M2 and N in Eq.(23) for random networks with Poisson distribution 
and networks with the normal distribution in degree distribution. From the figures, we observe 
a relation M2 ~ const., independently of N. The linear lines show N in the right hand side of 
Eq.(23). Though one intersection point in the Fig. 9 and Fig. 10 only makes Milgram Condition, 
M2 > N means that two degrees of separation is realized. It is getting quite difficult that the 
Milgram Condition is satisfied for small (fc). Though N < 1000 in our simulations, we can also 
speculate where is the intersection point between two lines for large (k) due to M2 ~ const, for 
Normal distribution. Moreover we notice the constant values are almost proportional to (fc)^. It 
is natural that such situation occurs because loop structures can be neglected in two degrees of 
separation. Thus we can approximately estimate the critical point that the Milgram Condition is 
satisfied for arbitrary (k). In this connection, we obtain (k) (a few x 10^) for ~ (a few x 10*^) 
which is about the population of USA, in a random network. We think that the value of (fc) 
is as large as it is unrealistic, because some social researchers estimated the average number of 
acquaintances of a person is 290 [27] .[281. [29] . 



Figure 9: M2 with (fc) = 10 and (k) = 20 for Poisson distribution 




Figure 10: M2 with ((fc),cr) = (10.5), (25.5), (25, 12), (50.12), (50,24) for Normal distribution 
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Figure 11: M2 and for 7 = 2, 3, 4 in SF networks 




Sr'F>at'S"ti on ISIi_jmt>er 



Figure 12: M2/N and Connection Ratio v.s. degree number of separation for SF network with 

7 = 2 



The behavior of M2 is rather different from the Normal and Poisson random networks for scale 
free networks with degree distribution P{k) ^ k^^ . (k) depends on the index 7 in the scale free 
networks produced by the configuration model. As N increases, so M2 increases quickly than TV 
as is shown in Fig.ll. The increasing rate is larger for the smaller index . Aoyama[2] has pointed 
out that critical index is 7 = 2 in two degrees of separation. Fig..ll, however, shows that 7 < 2 
can not realize the two degrees of separation. The result is different from Aoyama's one. 

For 7 = 2, Mp/N and the ratio of the zero components in in our formalism are given in 
Fig. 12 where p means the p—th degrees of separation ("separation number"). Evaluating the ratio 
of the zero components, 7 < 2 can not realize two degrees of separation but rather four degrees 
of separation is realized where 75 percent of all the nodes are connected. From (23), the critical 
velue of \ogiQM2/N would be an upper areas of the line a little smaller than zero. The Mp/N 
satisfying the Milgram condition in Fig. 12 also support this results. So, our formalism is thought 
to be available in order to analyze degree numbers of separation. We shall devote discussion on 
six degree of separation based on the formulation in a subsequent article pp]. 

5 Summary 

In this article, we provided a reformulation of the string formulation proposed by [2] to analyze 
networks. Fusing adjacent matrix into the formalism, we reformulate the string formalism. Ac- 
cording to it, we introduced a series of generalized q-th clustering coefficients. Their function is 
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not yet considered in this article and left as a theme for research in the future. Then we introduce 
the R matrix in the formalism developed in this article instead of A. The power of R plays central 
role in the analysis of this article. Every term of the expansion of i?" can be also interpreted 
graphically and it would make a projection for the future in the estimation of i?" that has complex 
expressions for large n. 

On the latter half of this article, we apply the formulation to some subjects in order to mainly 
check consistency with former studies. We first evaluated the clustering coefficient for typical 
networks studied well earlier. We could confirm a validity of our formulation by these in some 
degree. Lastly we applied the formulation to the subject of two degrees of separation. We find 
that the result is a little different from Aoyama's one [2]; the separation number is not two but four 
at 7 = 2. It is noticed that by analyzing the number of zero-components in A^, our results are 
rather supported. 

The following problems are yet left in future: 

1. Finding explicit expressions of i?" for arbitrary n by applying our formalism, especially 
diagrammatic construction. Then finding a general formula for arbitrary n from the series 
of expression. 

2. Revealing relations between p-t/i degrees of separation and A^, (fc) or < fc" >. More definitely, 
discovering the relations between p and N, (fc) or < fc" >. 

3. Revealing relations between p-th degrees separation and various loop structures, especially 

C(q). 

The last one will be discussed in the subsequent paper [20]. This article give a sort of general 
formalism to investigate above problems, including preliminary studies of them. 
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A i?6 g^nd Tri?" 

In this appendix we give an explicit expression and the expressions of Tr of R" for n = 1 ~ 6. 
R^ is obtained after straightforward but long tedious calculation. We divide it into the following 
four parts to brighten the prospects of the caluculation. 

— ^ Q^ijQ'jkQ'klQ'lmO'mnO'nf X A^fc Aj; Afcm A/^A^/ A^; Aj^ Afc„ A// A^^ A^^ Afc/ 

- ^ O^ifafkakiaimClmnClnf X A^fcA/; Afc^A^^A^/Aj^ AfenA^^ 

k,l,m,n 

— dijajkakiaimG'rniCiif X \k^jl^km^li^mf^jm^lf^kf 
j,k,l^m 

+ ^if^f^^klO'lrnO'rni^ik^fl^km^li^mf, 

k,l,m 

=R'[l],f + R'[2],f + R'[3]if + R'[%f, (25) 

where A^fe = 1 — Sik- Furthemore we divide R^[l]if into the following four parts to brighten the 
prospects of the caluculation. 

G-ijG'jk^kl^lm^mn^nf X ^ik^jl^km^ln^mf 

j,k,l,m,n 

= E (^ijC-ikO-kiaimamnanf X Ajfe Aj; Afe^ A(„ A^/ A,; Aj^ Afe„ A;/ Aj„ 

j^k.l,ra.n 

— aijajkakiaiidindnf X \k^jl^ln\f^kn^lf^jn^kf 

j,k,l,n 

~ E ^ij^jf^fl^lrnamnanf X Aj/ Aj; A/^ Ai„Aj; Aj^ Aj„ 
+ E '^ij^jf^'fl^'li^inO'nf^f^jl^ln^jn^mf, 

=i?'5[l, + 2],j + 3]i/ + R%4]if. (26) 
The four terms are respectively expressed as follows; 
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1],^ =[A%f + [A%f{4 - {h + kf)) + [AGA],f ih + kp) - {AG A, A^},f - [A^GA%f 
+ 2[A(G2 - 3G)Ay + 3j2(ivajkakf[A^]jk - + 

i 

-[^4],, - 2{[A% + [A^]f,) + [AGA],, + {{[A%f + {[A^]},f)) 

([A%f{{h - - 1) + 1 - - {[A%fY + X][A2]^,[A2],/(a,,- + a,/) 



|^[^^]//(2fc/ + fci - 5) + [A3]ii(2/ei + /e/ - 5) + - 3/ci - 3/c/) 

i?6[l,2]i/+i?6[l,3]i/ = -Ai/ + [A%f)+A[AGA\if - {A\G^ - 3G}i/ - {AGA,A%f 

-A[A%f - E a,ja,f(^{[A%ff + i[A%fr) + 2{[A%, + [A']j,) - {[A%, + [A']f,)^ 

+a,f (^-2[A%f[A%f + 2[A%f{ki + /c/ - 3) + 2j2<^vajf{[A%j + [A^jfj)^ , 

R^[l,4]if =[^=^],/A^/(([^=^],/)^ - 3[A%f + 2). (27) 
R^[2]if, R^[3]if and R^[4]if are respectively given by the following expressions; 

R'[2y+R^[3],f = a,f (^2[A%f - {{[A%, + [A']ff) 7{[A%, + [A^]ff) + 22[A%j 
+ ^[A%f[A%f + H[^%iki + [A%fkf) + Y,[A% {ajf + aij) 

3 

- 4Eaya,/([A2],,. + [A%^) - &{A\G}if - 2[AGA]if + {A,AGA}u + {A,AGA}fX 

3 / 

^'[4]i/ =aif [[A%f - [AGA\^f - {A\ G] + b[A%f - {[A%, + [A%f)^ . (28) 

By unifying all the terms, we obtain the full expression of R^ . It is too long and complex that 
we do not describe it here. Lastly we give the expressions of Tr i?" appearing in Eq. (7). 
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Tr(i?2) = 0, 
Tr{R^) = Tr{A^), 

Tr{R^) = Tr{A^) - 3Tr{GA'^), +2Tr{A'^) + TriC^ - G), 

Tr{R^) = Tr{A^) - 3Tr{GA^) + 6Tr{A^) - diag{A^)Tr{A^) + Ndiag{2A^G - A^), 
Tr{R^) = Tr{A^) + 6Tr{A^) - 5Tr{GA^) - iTr{A^) + Tr{A'^G'^) - 6Tr{A'^G) + 4Tr{A^) 

+ 2TriAGAG) - ^{auf - Y^[A%jaij[A%j +6j2aij[A%j + Yl aijajkaki[A%k 
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